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Abstract: A detailed study of half-BPS line operators of higher rank 4d A/" = 2 theory 

engineered from six dimensional Aj^-i (2, 0) theory on a bordered Riemann surface with 

' f— i ' full marked points is performed. Geometrically, each 4d UV line operator is represented by 

'^ an irreducible bipartite web formed by three junctions on Riemann surface, and such web 

^~< structure is called higher lamination. Algebraically, the space of UV line operators is iden- 

,^ tified with the integral tropical a coordinates of the corresponding PGL(N,C) local system, 

and the space of IR line operator is identified with the cluster X coordinates of SL(N,C) 

'"^ local system. The expectation value of UV line operator at Coulomb branch parameter- 

(^ ized by X coordinates is calculated, and the result is a positive Laurent polynomial in X. 

0\ Using the expectation values, we calculate the operator product expansion (OPE) between 

fvq the line operators, which is then represented geometrically by higher rank Skein relations. 

■^ We also calculate the Poisson brackets of these line operators, and Frenchel-Nielson type 

^^ coordinates are constructed for Higher Teichmuller space, etc. 
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1. Introduction 

The line or loop operators of quantum field theory play a very important role in probing 
the dynamics of the theory: the expectation value of the Wilson or 't Hooft loop of 4d 
gauge theory can be used to determine the phases of theory; The loop equations provide 
an important clue to the discovery of dual string theory [1], and Wilson loops are very 
important observables in the study of AdS/CFT correspondence [2, 3]; The expectation 
value of certain special Wilson loop is related to the scattering amplitude in 4d A/" = 4 
super Yang- Mills theory [4, 5], etc. 

We are going to study half-BPS line operators of a large class of 4d AA = 2 quantum 
field theories engineered using 6d Ai\f-i (2, 0) theory on a Riemann surface T, with various 
defects [6, 7, 8]. A lot of studies have been done for Ai theory ^ [9, 10], and they have a 
lot of interesting applications: 

• They are important observables of 4d theory. 

• A basis of them gives a nice maximal set of commuting Hamiltonian of the Seiberg- 
Witten (Hitchin) integrable system, which is very useful for the quantization [11, 12]. 

• They form a very interesting (quantum) algebra, and the Nekrasov instanton partition 
function is related to the quantization of this algebra [13, 14, 15, 16]. 

• The framed BPS states can be defined and calculated [10]. 

The 4d line operators of Ai theory are derived from wrapping half-BPS surface opera- 
tors of 6d theory^ on one cycle of Riemann surface, and the classification should be related 
to the classification of cycles on S. This picture is confirmed for Lagrangian Ai theory 
in [9], which shows that 4d line operator is classified by non-intersecting closed curves on 
S. Interestingly, such set of closed curves have been studied in great detail in mathemati- 
cal literature and called laminations [20]. Using lamination language, one can extend the 
classification to non-Lagrangian Ai theories [10]. 

This paper serves as an attempt to do a systematic study about line operators of 
higher rank theory defined using Riemann surface S with full defects (regular and irregu- 
lar). Naively, the generalization is straightforward: one wrap various 6d half-BPS surface 
operators on closed curves to get 4d line operators. This picture seems perfectly right if we 
only consider the usual Wilson lines built from weakly coupled gauge group, since there is 
indeed 6d surface operator labeled by arbitrary irreducible representation of SU(A^) gauge 
group [21, 22, 23]. So the classification using 6d picture agrees with the result using field 
theory analysis [24]. 

However, this can not be the whole story for the following two reasons: firstly, unlike Ai 
theory, there is strongly coupled matter system which carries non-trivial Coulomb branch 
operators, and there should be other non-trivial line operators charged on them as the 

^By An-1 theory we mean the 4d theories engineered using six dimensional An-i (2,0) theory. 
^Six dimensional An~i (2,0) theory are the low energy effective theory of N M5 brane, and the basic 
surface operator of (2,0) theory comes from M2 brane ending on M5 branes [17, 18, 19]. 



usual Wilson loop is not suitable; Secondly, it is shown in [25] that webs on S can also 
appear as the 4d line operators using 2d-4d correspondence [26]. 

We will confirm in this paper that the webs are necessary for the classification of 4d line 
operators. The important tool is the operator product expansion (OPE) between ordinary 
Wilson and 't Hooft line operator, since new line operators must appear in the OPE if they 
exist. Remarkably, those new line operators can be represented by webs formed by simple 
three junctions, and a complete classification of these three junctions will be given. Using 
these webs and the cluster algebra tools, we are going to make following main conjectures^: 

• The space of 4d UV line operators is identified with the integral tropical a coordi- 
nates of the moduli space of decorated PGL(A'', C) local system on S; The cluster 
X coordinates ^ of the framed SL(A'^, C) local system are identified with the IR line 
operators, which also parameterize the Coulomb branch. This is a local description 
as there are a lot of cluster coordinates systems related by quiver mutations. 

• The expectation value 1{C) of a UV line operator C on Coulomb branch is calculated 
and expressed as a positive Laurent polynomial in cluster X coordinate. The leading 
order term for a UV line operator with coordinates (oi, . . . , a„) is 

iA(/:) = n^r + ---- (1-1) 

Using the canonical map, one can define OPE between line operators: 

l{Ci)*l{C2) = J2clcAJ^), (1.2) 

c 

r r I r 

where the sum is finite and G^^j^^ is a positive integer, and C^.^^ ^ = 1. 

• We find higher rank Skein relations which give the OPEs in a geometric way. 

• Using the Skein relations, any 4d half-BPS line operator can be represented by an 
irreducible bipartite web on S, and those web structures are called higher laminations. 
This is a global description, and the space of higher laminations is in one-to-one 
correspondence with the integral tropical a coordinates described in first statement. 

Let me explain in some detail what each statement means. As we discussed earlier, the 
space of UV line operators in the A\ theory is identified with the lamination space. On the 
other hand, it is shown in [20] that the lamination space is in one to one correspondence 
with the tropical o coordinates of PGL(2, C) local system defined on S. The cluster X 
coordinate is related to the one cycle on Seiberg-Witten curve, which is then naturally 
identified with the IR line operator [10]. Remarkably, the cluster coordinates (a quiver) for 
the higher rank theory defined using full defects are found in a seminal paper by Fock and 



^These conjectures are confirmed for some explicit examples. 

*Tlie most important thing about the definition of cluster coordinate is a quiver, which is actually the 
BPS quiver. The tropical a coordinate is just a set of discrete numbers defined on each quiver node, while 
the cluster X coordinate is a set of complex number defined on each node. 



Goncharov [27], i.e. the quiver for A2 theory on once punctured torus is shown in figure. 1. 
Based on the analogy with Ai theory, it is straightforward to identify the space of UV 
hne operator as the tropical a coordinate of the corresponding quiver. Here we make the 
definition much more precise: they actually take values k/N with k integer and also satisfy 
some special constraints (see (3.26) and (3.27) for the detailed constraints of A2 theory.). 
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Figure 1: Top left: the BPS quiver for Ai theory on once punctured torus, here the quiver nodes 
on the edges labeled by same red lines are identified. There is a cluster X variable on each quiver 
node. Top right: the tropical a coordinate of a line operator L. Bottom: the OPE between Wilson 
and 't Hooft line, and a bipartite web must appear in the OPE. 

The second statement is established by studying line operator represented by closed 
curve on S: one can calculate its monodromy of SL(A^, C) local system using very simple 
combinatorial method [27], and the answer is always a positive Laurent polynomial^ in 
cluster X coordinates, which is the expectation value of UV line operator as a direct 
generalization of what is proposed for A\ case [10] (we also call it canonical map following 
[27].). For the closed curve L shown in figure. 1, the answer is 
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These canonical maps have two remarkable applications: firstly the tropical a coordinates 
of this line operator can be read from the exponent of leading order term, see figure. 1; 
Secondly it can be used to define the OPE between line operators: we simply multiply 
canonical maps of two line operators and the answer should be expanded into the sum 
of canonical map of other line operators. Using the OPE, we find that the usual Wilson 
loops are classified by the irreducible representation of the SU(A'') group, and the OPE 
between them has the same support as the corresponding tensor product of irreducible 
representations. This strongly suggests that our results are correct. 



^We call them Laurent polynomial even though the exponents are fractional. 



Most interestingly, by doing OPE between the usual Wilson and 't Hooft line operators, 
we find extra line operator which is naturally described as a web formed by two three 
junctions, see figure. 1. We classify all the three junctions of any higher rank theory, and 
any 4d half-BPS line operator can be represented by the webs formed by three junctions. 
These web structures are called higher laminations in direct analogy with the Ai case; 
Some webs on torus are shown in figure. 2. This description is independent of the cluster 
coordinates system and is the geometric picture for 4d line operators of higher rank theory. 
In general, the OPE between any line operators can be found geometrically using local 
Skein relations, see figure. 2 for A2 theory. Such Skein relations are also important for 
the classification of webs: we only consider irreducible webs which can not be reduced by 
Skein relations. 




Figure 2: A: Skein relations for A2 theory. B. Bipartite webs on once punctured torus of A2 
theory, and edges of three junctions with same label are connected. 



There are some compelling evidences supporting our results: Poisson brackets for the 
line operators of A2 theory on once punctured torus are calculated, and the answer is 
in perfect agreement with the underlying integrable system: there are three commuting 
Hamiltonians, etc. There are theories which can be engineered either using Ai or A2 
theory, and the canonical map of same line operator agrees with each other perfectly! 

This paper is organized as follows: in section II, we review some basic facts of line 
operators of Ai theory, and we emphasize a coherent picture of geometric and algebraic 
descriptions. In section III, we carry out a detailed study of line operators for theories 
defined on once punctured torus, and find the basic three junctions. Skein relations, etc. 
Section IV gives a brief description for line operators of theories defined on general Rieinann 
surface. Finally, a conclusion is given in section V. Technical details about how to calculate 
monodromy is given in appendix A, and we collect some useful canonical maps in appendix 
B. 



2. Ai theory: line operator and lamination space 

This section is a review about the hne operators for Ai theory. The geometric picture for 
the classification is the so-called lamination, and the algebraic picture is the tropical cluster 
algebra [20]. Physically, the geometric picture for the Lagrangian Ai theory is discovered 
in [9] using field theory , and then extended to the non-Lagrangian theory using lamination 
in [10]. We are also going to discuss the canonical map, OPE, etc, and all of which will be 
generalized to the higher rank case in the following sections. 

Let's first say a couple of words about the four dimensional M = 2 theory we have: 
we start with a six dimensional Ai (2, 0) theory and compactify it on a Riemann surface 
Mg^f,-^p., here bi denotes ith irregular singularity and pj denotes the jth regular singularity 
^. The IR limit is a four dimensional UV complete theory (conformal or asymptotical free 
theory). The UV deformations are encoded in the geometric description: the dimensionless 
gauge coupling is identified with the complex structure moduli of the Riemann surface; 
the mass parameters are the coefficients of the first order pole of regular and irregular 
singularity; the dimensional couplings are the coefficients of the higher order pole of the 
irregular singularity. The IR description can also be easily found from this geometric 
description: a Hitchin equation with specific boundary conditions at the singularity is 
defined on Mg^f,i,pj, and the spectral curve of the Hitchin moduli space is identified with 
the Seiberg-Witten curve of 4d theory. 

The weakly coupled gauge group description '^ is described by degeneration limit of 
the Riemann surface: the long tube in the degeneration limit is associated with the gauge 
group and other components are the matter system . The basic building blocks (matter 
contents) are: 

• Tri-fundamental represented by a sphere with three regular singularities. 

• Dtv type® Argyres-Douglas theory represented by an sphere with one irregular sin- 
gularity and one regular singularity . 

The first piece is free while the second one is a strongly coupled isolated SOFT, see figure. 
3A. Each regular singularity carries a SU(2) flavor symmetry, and the weakly coupled gauge 
group description of the four dimensional theory is derived by gauging these SU(2) flavor 
symmetries, see flgure. 3B. Different duality frames correspond to different decomposi- 
tions of the full Riemann surface into the above two pieces, which is a remarkably simple 
description of S duality [6] . 

For various applications below, we need to replace an irregular singularity with a 
boundary plus marked points whose number depends on the type of singularity, which 
can be seen by studying the Stokes phenomenon around the irregular singularity, see [28]. 
Therefore, we are going to consider a bordered Riemann surface Mg^{|.^„, here g is the genus 



®The irregular singularity means that the Higgs field of the Hitchin equation has higher order pole, while 
the regular singularity means first order pole. 

^For the sphere with one irregular singularity or one irregular and one regular singularity, one has 
strongly coupled Argyres-Douglas theory. 

*It is called Dat type theory since Dn quiver appears in the mutation class of its BPS quiver. 
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Figure 3: A: The two basic matter systems for four dimensional theory built from six dimensional 
Ai theory: the tri-fundamental represented by a sphere with three regular singularity and the D 
type Argyres-Douglas theory represented by a sphere with one irregular and one regular singularity. 
B: In one weakly coupled gauge group duality frame, the theory is formed by gluing above two pieces 
together: physically this is achieved by gauging the diagonal flavor symmetries. C: We replace each 
irregular singularity with a boundary with marked points. 



of the Riemann surface, bi is the number of marked points on the ith boundary, n is the 
number of bulk marked points, see figure. 3C. The Coulomb branch dimension of the 4d 
theory is 



rir 



n + 






+ 35-3, 



(2.1) 



\b^-■i^ 



where [^^] means taking the integer part of the number inside the bracket. 

Each regular puncture contributes a mass parameter (as it carries a SU{2) flavor 
symmetry.), and each boundary with even number of marked points also contributes a 
mass parameter (since it carries a U{1) flavor symmetry.), and the rank of the charge 
lattice is 

i? = 2n,. + n/ = 3n + ^ 6i + 65 - 6. (2.2) 

i 

We also pick an anti-clockwise orientation for each boundary (the bulk puncture can be 
thought of as a boundary without any marked point.). A bordered Riemann surface with 
anti-clockwise orientation for each boundary is the geometric avatar for our study of line 
operators. 



2.1 UV and IR line operator: geometric picture 

The 4d theory always has a Lagrangian description for bordered Riemann surface M^^o.rn 
so one can use field theory tools to study its half-BPS line operator. For a four dimensional 
M = 2 gauge theory with gauge group G, the Wilson operators are labeled by irreducible 
representations R of G and take the form 

Wr = Trj^P / exp(iA + Recjyds), (2.3) 

where i;^ is a complex scalar in the vector multiplet and ds is the line element determined by 
the metric. If the integral is taken along a straight line or a circle ^ in flat -R4, this Wilson 
line preserves a half of AA = 2 supersymmetry. Equivalently, the Wilson line operator is 
classified by the highest weight in the weight lattice Ac of the gauge group G. The 't Hooft 
operator is a disorder operator and defined by specifying the singular behavior of various 
fields around the loop or line in the space time: 

F = smedeAd(j), (p^i — , (2.4) 

2 2r 

and the line operator is classified by the highest weight in the coweight lattice in A^^,, or 
the weight lattice of Langlands dual group G. For general Wilson- 't Hooft line operators, 
Kapustin found that the line operators are classified by the following data [24]: 

(^, u) G kcw X A^, (/i, u) ~ (w./u, Lo.v), uj £W, (2.5) 

where W is the Weyl group. The allowed magnetic weight has to satisfy the Dirac quan- 
tization condition in the presence of matter fields. Moreover, when the flavor symmetry 
is available, one can introduce the background gauge field for these flavor symmetries and 
consider the line operator carrying magnetic charge on it. 

Let's apply the above general classification to the theory defined by Riemann surface 
Mgfi^n, namely a Riemann surface with n regular punctures. There are 3g — 3 + n gauge 
groups and n mass parameters, and a duality frame is described by a pants decomposition 
of the Riemann surface (since there is only tri-fundamental matter.). The line operator is 
labeled by {6g — 6 + 3n) integers 

{pi,P2, ■ ■ ■ ,P3g-3+2n; Ql, 92, • • • , QSg-S+n)- (2.6) 

The Weyl transformation acts as {pj, qj) — )• {—pj, —Qj) for j = 1, 2, . . . , 3g — 3 + n, and 
Pj — )• —pj for j = 3g — 2 + n, . . . ,3g — 3 + 2n. Using Weyl transformation, we may take 
p to be positive, and they are the magnetic charges, so they have to satisfy the condition 
Pi + Pj + Pk ^ 2Z for three gauge groups coupled to a single pants, since the matter from 
the pants has to satisfy the Dirac quantization condition. 

Using above field theory results, it is shown in [9] that the space of line operators are 
identified with the space of self- and mutually nonintersecting unoriented closed curves 



^We call these two types as line operators for simplicity, as the space-time shape is not important for 
our later purpose. 




Electric charge: Winding number 



Magnetic charge: intersecting numbers with boundary 



Figure 4: The Dchn-Thurston coordinates of a set of closed curves are defined using a pants 
decomposition: the pi coordinates are defined by simply counting the weighted intersecting number, 
while the qi coordinates are defined by counting the oriented winding number. 

on Mgfl^n- Therefore, instead of studying the line operator from 4d field theory point of 
view, we can just study these closed curves on the Riemann surface. Now let's start with 
a collection of closed curves (non-intersecting) with positive integer weights ^^ , and one 
can construct a set of Dehn-Thurston coordinates based on the pants decomposition: the 
coordinates pi (always positive) on a boundary circle of the pants is defined by counting 
the weighted intersection number, and the qi coordinates (it can be positive and negative 
depending on the winding direction) on a tube are defined by counting the winding number 
of the curves, see figure. 4. On the other hand, given a set of coordinates based on pants 
decomposition, one can easily find a set of closed curves by doing the inverse construction: 
one put Pi marked point on the i th boundary of the pants, and connect them pairwise 
without intersection on each pant (which is possible if they satisfy the condition pj+pj+p^ G 
2Z on each pant, and the construction is unique), and also the number of winding on the 
tube is determined by the coordinates qi , and we construct a set of non-intersecting closed 
curves with weight. Finally we combine the closed curves in same homotopy class, which 
will give back a geometric objects. 

The above Dehn-Thurston coordinates are directly related to the field theory descrip- 
tion [9], however, it has several disadvantages: firstly, the transformation of the coordinates 
is quite complicated under the change of pants decomposition corresponding to the S dual- 
ity; Secondly, there is no pants decomposition for the Riemann surface with boundaries and 
the above definition of the line operator is not applicable. In next subsection, we are going 
to introduce another set of coordinates based on the triangulation of the Riemann surface, 
and the coordinates has simple transformation rule under the change of the triangulation. 



^"The label is taken as the charge of that line operator in the duahty frame where the closed curve 
represents a weakly coupled gauge group. 



moreover, the definition can be extended to all bordered Riemann surface. 

2.1.1 UV line operators: lamination and its coordinates 

The set of the non-intersecting closed curves are called laminations in the mathematics lit- 
erature [20] , so space of line operator of 4d theory is naturally identified with the lamination 
space. Let's now introduce the precise definition of Lamination on bordered Riemann sur- 
face. Before that, let's give the definition of special curves: a curve is called special if it 
can be contracted to a bulk marked point or a boundary marked point, and a curve is 
contractible if it can shrink to single point which is not a marked point. Now let's spell 
the definition of the Lamination: 

Definition: A integral A-lamination on a bordered Riemann surface is a homotopy 
class of a collection of finite number of self- and mutually nonintersecting unoriented 
curves either closed or connecting two points of the boundary disjoint from marked points 
with integral weights and subject to the following conditions and equivalence relations. 

(1) Weights of all curves are positive, unless a curve is special. 

(2) A lamination containing a curve of weight zero is considered to be equivalent to 
the lamination with this curve removed. 

(3) A lamination containing a contractible curve is considered to be equivalent to the 
lamination with this curve removed. 

(4) A lamination containing two homotopy equivalent curves with weights u and v is 
equivalent to the lamination with one of these curves removed and with the weight u + v 
on the other. 







a'=l 



Figure 5: A: Laminations on an annulus with one marked points on each boundary, which represent 
the pure SU(2) theory. B: Constructing coordinates for the lamination using the triangulation. 

Notice that the weights for the curves around the marked points can be either positive 
or negative, while the other curves always have positive weights. An example of lami- 



10 



nation on an annulus with one point on each boundary is shown in figure. 5. Instead 
of constructing coordinates using the pants decomposition, we use the ideal triangulation 
whose vertices are the marked points. An ideal triangulation is formed by simple curves 
whose ending points are the marked points, equivalently one can use a quiver to encode 
the triangulations of the internal edges, and this quiver is actually the BPS quiver [29]. An 
important fact is that the number of internal edges in the triangulation is 6g — 6 + b + 3p 
which is exactly the same number of the charge lattice We give some example of triangula- 
tions and quivers for various bordered Riemann surface in figure. 6, more details about the 
ideal triangulation can be found in [30]. Given a lamination and a triangulation, we can 
construct the coordinates (labeling) for the lamination using the following simple formula: 



ttE 



IE 



UJi 



(2.7) 



where the sum is over all the curves intersecting with the curves, see figure. 5B. The 
coordinates are called tropical a coordinates and taking half-integer values. 






Figure 6: Top: The ideal triangulations of bordered Riemann surface: disc with six marked points, 
sphere with four punctures, and annulus with one marked points on each boundary. Bottom: the 
corresponding quiver. 

On the other hand, given a set of constrained half-integers (the constraints will be 
spelled later) on each edge of the triangulation, we can construct a lamination. First of 
all, we can assume all the coordinates are positive, since if some of the coordinates of aj is 
negative, we can first shift all the coordinates by a positive integers u such that a^ = Ui + oj 
are all positive. If we can construct lamination L for a^, we can add closed curves with 
weight —uj around each marked point to L (these curves will not intersect with L ) to get 
the lamination L for the coordinates Oj. So let's start with positive half-integer coordinates 
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ue which satisfy the following two conditions for the coordinates on a triangle: 

I «a — 0/3 1 < a-y < CLa + a^ 

da + ai3 + a^ G z. 



(2i 



Now for each edge in the triangle < ABC >, we put 20^ number of points on it, and 
for a triangle, we connect these points pairwise such that there is no intersections. The 
construction is unique, since we can divide the marked points on each edge into the two 
parts as (and the coordinates on each edge are {AB, BC, CA) = [d^, dp^da))'- 



a + b = 2a^ 
a + c = 2da 
6 + c = 2a^ 



a = da + d^ — di3 
b = di3 + d-y — da 
c = da + dp — d^, 



(2.9) 



and the solutions are ^^ positive integers due to the constraint (2.8), and they determine 
how to connect the points on the boundary of the triangles: i.e. there are a curves around 
vertex A. The ideal triangulations are not unique, and different triangulations are related 




Figure 7: A: Given a set of tropical a coordinates on a triangle, one can construct a set of non- 
intersecting curves, here the coordinates are {AB^ BC, CA) = (5, 3, 4). B: The rule for determining 
the tropical x coordinate. 

by the local moves called flip (see figure. 5). The a coordinates change in a very simple 
way under the flip: 

Cg = — ttg + max{aa + flci o-b + o^). (2.10) 

If there is a boundary, one need to put further constraints on the a coordinates: the 
coordinates on the boundary edge of the triangulation should be zero. The reason for this 
further constraint is that the external edge is not included in our description of A/" = 2 
theory, as the number of internal edges are exactly equal to the rank of the charge lattice. 



^^The number of curves around a vertex is equal to the sum of the points on two edges connecting this 
vertex minus the number on the opposite edge. 
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One might be curious about the strange constraint (2.8), and the reason is that these 
coordinates are actually the tropical a coordinate of the PGL(2, C) local system, which is 
the Langlands dual group of SL(2, C). We assume that there is a SL(2, C) Hitchin system 
on the Riemann surface for the 4d field theory description ^^ , and the line operator is 
taking value on the tropical a coordinates of PGL(2, C) local system. For more discussion 
on the global issue about the gauge group, see the appendix of [7]. 

One can construct another set of coordinates for the lamination ^^. The definition is 
the following: pick a vertex of the diagonal edge of the quadrilateral and name it as A, 
then label the other vertices as B,C,D in clockwise direction; if a curve goes through AB 
and CD, its coordinate is positive, and the coordinates is negative if it goes through AD 
and BC, see figure. 7. One can check that the x and a coordinates are related by the 
following simple formula 

Xi = -eijOj. (2.11) 

where Cij is the antisymmetric tensor read from the quiver. These x coordinate is always 
integer, and can be both positive and negative even for the non-trivial closed curves. 

2.1.2 IR line operator: one cycle on Seiberg-Witten curve 

The UV theory is described by compactifying 6d theory on a punctured Riemann surface, 
and the lamination on Riemann surface T, will give us the UV line operator. From M 
theory picture, we have a 4d theory on two M5 branes, and M2 brane ending on M5 branes 
(wrap on lamination on S) gives 4d line operators. What about the IR line operator? In 
the IR limit, these multiple M5 branes become a single M5 brane compatified on another 
Riemann surface 'E,^ which is identified with the Seiberg-Witten curve [31]. It is natural to 
expect that M2 brane wrapped on various one cycle of S'^ will give us IR line operators. 

There is actually a distinguished basis on these one cycles based on the ideal trian- 
gulation of the Riemann surface S. The construction goes as follows: first of all, we can 
replace the triangulation by a bipartite network as suggested in [32] . Usually one associate 
a face variable Xi to each face in the original bipartite network, and one can read a quiver 
from the network. Given a bipartite network on a Riemann surface, one can construct 
a conjugate Riemann surface T,'^ using the so-called untwisting procedure, which is quite 
familiar in the dimer study, see [33]. The details are not important for us, and we want to 
point out the crucial point: since the boundary of the face of network on S becomes the 
one cycle on T,^, which will give a IR line operator. It is then natural to interpret the face 
variable X as the expectation value of the IR line operators . On the other hand, these 
one cycles on the SW curve is also used to parameterize the Coulomb branch through the 

^^If the gauge group of the Hitchin system is PGL(2, C), then the hne operators will take values in the 
tropical a coordinates of the SL(2, C) local system and will take integer values. 

^^This is called the X lamination, which is also described by the curves on Riemann surface, but the 
precise definition is quite different from the A lamination defined earlier. The crucial point is: every A 
lamination is a also a X lamination, so one can define its tropical x coordinate [20] . The contrary is not true 
though: there are X lamination which is not A lamination, namely one can include the curves connecting 
the punctures. 

^^Classically, the expectation value of the IR line operator is the length of that one cycle. 
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famous Seiberg-Witten construction, so the cluster X coordinates not only represent the 
IR line operator but parameterize the Coulomb branch of 4d theory! 

2.2 UV and IR line operator: algebraic picture 

Geometrically, the 4d line operator is related to the weighted closed curves on the Riemann 
surface. Its tropical a and x coordinates can be constructed using the triangulation and 
their transformation property under the change of triangulation is simple . In fact, the 
above picture can be put into a general algebraic picture called (tropical) cluster algebra 
[34, 20]. 

The (tropical) cluster algebra consists a family of seeds (ejj,Xj, Oj), here e^j is an 
antisymmetric matrix which can be represented by a quiver, and Xj and Oj are rational, 
real or integral numbers ^^ associated with each quiver node. Different seeds are related 
by the so-called mutations. The mutation acting on the quiver is 



-eij \i i = k or j = k 

"'^ ^ ejj +sgn(eifc)[eifcefcj]+ otherwise 



4= , 7"^. . ..-..,.-. (2.12) 



The quiver mutations can be represented beautifully using the quiver diagram: A quiver is 
a directed graph where multiple arrows between two vertices are allowed, which is derived 
using €ij as follows: attach a quiver node for i = 1, . . .n, and there are €ij arrows between 
node i and node j ^^. The quiver mutation for a quiver without one and two cycles (such 
quiver is called 2-acyclic) is defined as the following: Let Q be a quiver and k a vertex of 
Q. The mutation HkiQ) is the quiver obtained from Q as follows, see figure. 8: 

1) for each sub quiver z — )■ fc — )■ j, create a new arrow between ij starting from i; 

2) we reverse all arrows with source or target k; 

3) we remove the arrows in a maximal set of pairwise disjoint 2-cycles. 





t rs-t 

Figure 8: The quiver mutation. 

In changing the triangulation, the quiver is changed as quiver mutation (Seiberg du- 
ality), and the coordinates are changed in following way 

Ofc = max{ak[eik]+, ai[-eik]+) - au (2.13) 

-Xk if j = k 

Xj + efcj?nax(xfc,0) if ekj > (2-14) 

Xj + €kjmax{—Xk, 0) if ekj < 



Xj 



'^^The integral numbers does not necessarily take integer values, but they should be some discrete numbers. 
^®If eij is positive, the quiver arrows are pointing from node i to node j; otherwise, they are pointing 
from node j to node i. 
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The ordinary cluster algebra consists of also triples {€ij,Xi, Ai) which are related by quiver 
mutations, where eij is the same as the one used above and it has the same form as quiver 
mutations. The transformation of the X coordinates under mutation is 



' _ j ^k^ if J = k 



^k 

X,{l + X-'"'''^'")-'^>'ifj^k, 



and the A coordinates transform in following way 

^^ = _Ai ±11^ . (2.16) 

" Ak 

A two form can be defined on A space 

to = eijdlogAiAdlogAj; (2.17) 

and a degenerate Poisson structure can be defined on the X space: 

[Xj,Xj] = eijXiXj. (2-18) 

These structures are compatible with the cluster transformation, i.e. if you express {Xi, eij) 
in terms of (X^, e^ •) using the cluster transformation rule, and you will get the same form 
expressed in terms of {X^, e,). 

The above two types of coordinates are related by the so-called tropification procedure. 
Here let's describe a little bit about the idea of tropification, i.e. the triple (ejj,Xj,aj) can 
be thought of as the tropical limit of {eij,Xi, Ai). Let's explain a little bit of tropification: 
the tropification about a polynomial 3a^ + a^b + 1 is max(3a, 2a + 6, 0); more generally, for 
a rational function without subtraction, the tropification is 

TVop(/(ai, 02, ... , an)) = lim,^^€logif{e"'\e"'^,. . . , e^"")). (2.19) 

It is then easy to see that {€ij,Xi,ai) is a tropification of (€ij,Xi,Ai) by comparing their 
transformation property under the mutation. 

2.3 Canonical map: relating UV and IR line operators 

Classically it is natural to identify the expectation value of a UV line operator as the 
length of the closed curves on the Riemann surface, and the length should depend on 
the Coulomb branch, which is parameterized by the cluster X coordinates. There is a 
remarkably simple method to calculate the length: it is related to the monodromy around 
the curve in the SL(2,C) local system, and the monodromy can be calculated using very 
simple combinatorial methods [20]. Basically, one first construct a dual graph from the 
ideal triangulation: put a trivalent vertex inside each triangle and glue them together using 
the triangulation. Then put a square around each edge and a hexagon around each vertex, 
and this structure is called as monodromy graph^^. The orientation of edges in monodromy 



^^This prescription is different from what is usuaily written in the iiterature, and we adopt this definition 
because it is immediateiy generalized to higher rank. 
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graph is chosen as fohowing: the edges on the cycles around the marked points are taken 
to be in anti-clockwise direction, and the edges in hexagon is required to be in clockwise 
direction which can be used to fix the orientation of remaining edges. 

There are three types of edges in the monodromy graph: S edge which intersects with 
the dual graph, E edge which is parallel with the dual graph, and finally the F edge which 
is inside the triangle, see figure. 9. We associate the following matrices with the three 
types of edges: 



S 



-1 

1 



,E 



1 



.F 




If a path is going opposite direction to the prescribed orientation for an edge, we simply 











M 


c; 


y 




c 


^ 


r 






a 


L 





S: Cross edge of dual graph 



F: Inside triangle 







Figure 9: Top: the dual graph and the monodromy graph are built from the triangulations. 
Middle: One triangulation has three edges labeled by a, b, c, and the opposite edges are identified. 
We also draw various closed curves, and the quiver. Bottom: the Skein relations. 

use the inverse matrix. It can be checked that the monodromy around the square and the 
hexagon is an identity matrix up to a scalar, therefore it defines a PGL(2) local system. 
Now the monodromy of a closed curve C is found as following: find a closed path on the 
monodromy graph which is in the same homotopy class as C, and then multiplying the 
matrices of segments in this path which will give the monodromy. For our purpose, a SL(2) 
local system ^^ is required and the monodromy M is defined as 



M = |det(Mo)-^/2Mo| 



(2.20) 



here Mq is the monodromy using the above prescription and we take the absolute value to 
make all the expressions positive. Using the monodromy, the expectation value of the UV 
line operator or canonical map ^^ for curves in lamination are 



^^Since we assume the gauge group is SU(2). 

^^Since the expression map tropical a coordinates to an expression in cluster X coordinates. 
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For a non-special closed curve / with weight k, the canonical map is I(/) = Tr(M ). 

For a special curve / with positive (negative) weight k, the canonical map is to take the 
highest (lowest) eigenvalue of M' '. For example, the eigenvalues of the monodromy 
matrix have the following form (winding once): 



and we define the canonical map of curve with weight one as I(/^) = W X^ and the 

_ 1 
curve with weight minus one as I(/~) = H "'^j ^ • 

• For a curve connecting two points on the boundary, the canonical map is defined as: 

I(0 = (l0)*M*rj, (2.22) 

here the definition of M is subtle as there are two choices for the starting and ending 
edges of the path / in the monodromy graph. Our prescription is that the starting 
and ending edge of / has the same orientation as the curve. 

• For the line operator with coordinate, the canonical map is 1(0) = 2. 

• For a lamination / consists of many curves /j, and the canonical map is defined as 

l{l) = Wl{k). (2.23) 

i 

Example: Let's apply the above discussion to the theory defined by an once punctured 
torus, which represents the so-called SU(2) M = 2* theory. The triangulation and the 
quiver is shown in figure. 9, and there are three quiver nodes. There are three non-trivial 
homology cycles which we call L, M, C, and the line operator associated with this three 
cycles has tropical a coordinates (^, |, 0), (0, ^, |), (|, 0, |), and the canonical map is found 
to be 



I{L) = ■ 


1 b 
Vab V a 


^''''"^-Jii^fi^ 


IiC) = 


1 [^ 
\Jca V c 


/^, I(L + M) = ^^ + 
a 


I(P+) = 


--ahc, I{P-)-- 


= a-'b-'c-\ 



be 



ac , , — bJlic Jac 

— + b^/ac + -^ — + ^ 

ao a abc 

(2.24) 

here L + M is the cycle shown in figure. 9. One can calculate the OPE using the canonical 
map: let's multiple the canonical map of two line operators, and expand the answer in 
terms of the canonical map of other line operators. A simple example is 

I{L)*I{M) = I{L + M) + I{C), (2.25) 

Geometrically, the OPE can be found using the Skein relation, see figure. 9, and the above 
OPE can be easily found using Skein relations. Here are some properties for the canonical 
map I{1) (for more details, see [27]): 
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• The canonical map is a Laurent polynomial in Xi with positive integer coefficient, 
and the leading order term for a line operator with coordinates Oj is Y\ ^f ' • 

• There is an OPE between line operators I(/i)I(/2) = ^iCi i I(/), here Cf ; is a 
positive integer, with C^^i ^ = 1. 

• The Poisson structure between line operators can be easily calculated using the canon- 
ical map and the Poisson structure of cluster X coordinates. 

2.4 Summary 

Let's summarize the important facts about Ai line operators (see figure. 10): 



4d UV line operators 



4d IR line operators 



.Geometric: Closed curves on Riemann surface 



Equivalent 



■ Algebraic; Tropical a coordinate 



Canonical map 



— > Cluster X coordinate 



Figure 10: The basic picture of Ai line operators. 

Each 4d UV line operator is represented by a lamination, which is a global description. 

One can parameterize the UV line operator using the cluster algebra and tropical a 
coordinates: given a lamination, one can construct its coordinates; and given a set 
of coordinates, one can construct a lamination. This description is local as there are 
a lot of ideal triangulations. 

The space of IR line operator is identified with the cluster X coordinate. 

There is a canonical map 1(1) between the lamination space and the cluster X coor- 
dinates, and there is an OPE between them, which has very interesting property. 
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3. An-i theory: torus with one full puncture 

We are going to generalize the story of Ai line operators to ^at-i theory defined using 
full regular and irregular singularities^^. Similarly one can replace the irregular singularity 
with a boundary with full marked points, and the geometric object is still a bordered 
Riemann surface (see [28] for the exact correspondence between the bordered Riemann 
surface and the physical 4d A/" = 2 theory). In the Ai case without the boundary, the field 
theory tools are available since the theory has Lagrangian description, which provides a 
lot of informations. Now although the weakly coupled gauge group description is still valid 
and can give us some clues, there are always strongly coupled matter and it is clear new 
phenomenon must appear. 

On the other hand, the geometric picture we learned from Ai case is still valid, and 
closed curves on Riemann surface should still give 4d line operators. However, it should be 
clear that this is not enough as each closed curve give a weakly coupled gauge group, and 
these line operators can not be used to probe the strongly coupled matter part. 

We do have the cluster algebra tools on hand, so it is straightforward to argue that the 
UV line operators are parameterized by the integral tropical a coordinates of the cluster 
algebra, and the IR line operator is given by the cluster X coordinate. However, the specific 
tropical a coordinates (what kind of constraints one need to impose on the coordinates) 
are not known, and given a set of tropical coordinates, we do not have a geometric picture. 

Given all these difficulties, it is probably useful to look at the simplest theory: torus 
with one full puncture. It turns out that we are going to discover the generic features about 
the line operators of Ajv_i theory, and remarkably we will establish the same story as the 
Ai theory. The key is the canonical map of the line operator associated with closed curve 
and the OPE! 

3.1 A2 theory 

3.1.1 Closed curves: ordinary Wilson line 

Let's consider 4d AA = 2 theory defined by compactifying six dimensions A2 theory on a 
torus with a full puncture. The field theory in one duality frame is derived by gauging the 
diagonal of two SU(3) flavor symmetries of T3 (Eg) theory. This theory has three Coulomb 
branch operators: two are carried by SU(3) gauge group, and one is carried by the T3 
theory, see figure. 11; there are two mass parameters carried by the full puncture, so the 
charge lattice has dimension 8. 

The cluster coordinates (or the quiver) for this theory is constructed as following: one 
start with a triangulation of once punctured torus, which is the same as the Ai case, and 
the new features are: there are 2 nodes on each edge, and one extra nodes inside each 
triangle. The quiver is shown in figure. 11. 

It is natural to identify the integral tropical a coordinates as the space of half-BPS 
line operators. Although this statement is quite powerful, we actually want to know the 



^"in the regular singularity case, the full puncture means the most generic singularity, and the specific 
irregular singularity which will give the full marked points is explained in [28]. 
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Figure 11: Top: The 4d theory defined by putting 6d A2 theory on a torus with one fuU puncture, 
and it is derived by gauging the flavor symmetry of T3 theory. Bottom Left: The cluster coordinates 
for A2 theory on once punctured torus; Bottom Right: the closed curve L is represented by the red 
path in the monodromy graph. 

geometric counterpart which is more useful in many contexts. Namely, we would like to 
find the geometric representation for each line operator labeled by a tropical a coordinate. 

We do have a natural geometric object on the Riemann surface: the closed curve. Since 
each closed curve represent a weakly coupled gauge group, it is natural to expect that one 
should get a Wilson loop by wrapping six dimensional surface operator on it, and we would 
like to know its tropical a coordinates. According to A\ theory, the tropical a coordinates 
can be read from the top component of the canonical map calculated from monodromy; 
We expect this fact is also true in the higher rank theory, and luckily the way of calculating 
the monodromy is given in [27]. We are going to use the monodromy calculation to find 
the tropical a coordinates for each closed curve. 

There are three simple nontrivial cycles which we denote as L, M, C which can be 
identified as three duality frames of the corresponding gauge theory, see figure. 9. Let's 
consider one duality frame in which the L cycle represents the SU(3) gauge group. A closed 
curve wrapping around the L cycle should give us a Wilson loop. The monodromy of the 
closed curves can be similarly found using the following three by three matrices in the same 
monodromy graph as used for A\ theory: 



S[3] 




E[o,3] 



aofli 
ai 
1, 



Face[x, 3] 




The above prescription actually gives a PGL(3) matrix and SL(3) canonical map is given 
by the following formula 

I(/) = (Det(M))^5Tr(M). (3.1) 
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The canonical map of L (we call it as Li.) is found using the path in monodromy graph 
shown in figure. 11; and the monodromy matrix is 

M = E[a, 3] * Face[2;, 3] * E[6, 3] * S[3]"^ * Face[y, 3]"^ * S[3]-^ (3.2) 

Using (3.1), the canonical map for Li is 

... 1 , ^r , -^'^^r , xV3aV35V3 

^ ^' 2/3 1/3 1/3 2/3,1/3,2/3 "^ 2/3 1/3 V3 2/3,1/3 ^ 1/3 1/3 2/3,1/3 ^ 1/3 l/3,l/3~'~ 

1/3,2/3,1/3 1/3 1/3,2/3,1/3 ^l/3„2/3^1/3,2/3,l/3 , , , , 

^ ' % Pi , ^ ' "■! % Pi , X ' y ' a^ bQ b^ 1/3, 2/3„2/3 1/3,2/3, 1/3 .„ „^ 

1/3 1/3 2/3+ 1/3 1/3 + 1/3 +^ y «0 «1 ''O ^1 • l-^-^J 

This is a positive Laurent polynomial with positive integer coefficient, and the exponents 
of the leading order term should be identified with the tropical a coordinates of Li^^: 
(ao,ai, bo, 61, cq, ci,x,y) = (|, |, |, ^, 0, 0, ^, |). We now argue that Li should be identified 
with the Wilson line in fundamental representation of SU(3) gauge group by looking at its 
dual tropical x coordinate: 

Xi = -ttjCij, (3.4) 

and the answer for Li is (—1,0,1,0,0,0,0,0). Now take one of the edge (say b edge) of 
the triangulation, and there are 2 dots on it which can be thought of as the Dynkin nodes 
for the A2 group, and our line operator has coordinates (1,0) which is exactly the Dynkin 
label for fundamental representation of SU(3) group, so it is natural to identify it as the 
Wilson line in fundamental representation, see figure. 12. 

The Wilson loop in anti-fundamental representation is naturally identified with the 
closed curve with opposite orientation (this is different from Ai case in which the orientation 
does not matter, since the fundamental and anti-fundamental representation of SU(2) group 
is the same.) The canonical map is calculated using the matrix M~^, and the answer is 

^^^'^ ^ xV3y2/3„2/3^1/3^2/3^1/3 + ^1/3^2/3^2/3^1/3^1/3 + ^l/3„2/3^'j/3^1/3 + ^1/3^1/351/3 + 

„1/3a1/3/,2/3 i/3„1/3,1/3,2/3 ^2/3„1/3„1/3, 1/3, 2/3 , , , , 

y "0 "l , y % "0 "l , X ' y ' a^ fep b^ 2/3, l/3„l/3 2/3, 1/3,2/3 .„ .s 

1/3 2/3 1/3 + 1/3 1/3 "^ 173 + ^ y «0 "1 "O Ol • l-^-^J 

Now a coordinates for this line operator is (ao, ai, bo,bi, co,ci,x, y) = (^, |, |, |,0, 0, |, |), 
and the dual x coordinates is (0,-1,0,1,0,0,0,0). The x coordinates on edge b is (0, 1) 
which is exactly the Dynkin label for anti-fundamental representation. 

The other irreducible representations for A2 group are labeled by two positive integers 
{m,n). In group representation theory, these other representations can be found by doing 
the tensor product of m fundamental and n anti-fundamental representation, and expand 
them in terms of irreducible representation in which {m, n) is the leading order component. 
Here we use exactly same strategy to find the line operator for the other representations: a 



^^We use the same letter say ao to denote the cluster X coordinates and the tropical a coordinates, as 
the tropical a coordinates can be read from the exponent of the monomial in cluster X coordinates. 
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Figure 12: Left: the tropical a coordinates for line operators in cycle L with two different orien- 
tations. Right: the tropical x coordinate for L\ and Li- 

line operator with Dynkin label (m, n) can be found as leading component in the OPE of 
product L™L2 , since the leading order component in the OPE has the tropical x coordinates 
(m, n) on edge h. 

The tropical a coordinates of the line operators in representation (m,n) is ma{Li) + 
na(L2), and the canonical map can be found by carefully doing the OPE. A very useful 
fact is that the canonical map for the line operator in (n, 0) or (0, n) representations can be 
found using the monodromy matrix M" and M~" respectively. To find the canonical map 
of the general representation {m,n), one need to do the OPE expansion of the product 
L™L2. Not quite surprisingly, the OPE is closely related to the tensor product of the 
corresponding representations. More precisely, the support of the the OPE is exactly the 
same as the corresponding tensor product expansion of the group representation theory, 
but the coefficient of our OPE is not equal to the multiplicities of the tensor product. Here 
are some of the examples we found: 



I(wi) * l{u}2) = I(wi + u;2) + 1(0), (3 «) 3 = 8 + 1) 

I(wi) * I(a;i) = I(2wi) + 21(^2), (3 «> 3 = 6 + 3) 

I(wi + U2) * I(wi) = I(2a;i + U2) + 2I(2a;2) + 2I(lji), (8 3 = 15 + 6 + 3). 



(3.6) 



Here we use wi to denote the line operator in fundamental representation (Li) and UJ2 to 
denote the anti- fundamental representation line operator (L2), and the trivial line operator 
(with zero a coordinates) give a number 3 ^^ . Let's now prove the second OPE in above 
formula. The canonical map for (2, 0) line operator can be found using the monodromy 
matrix M^, and the result is given in the appendix B. Using the canonical map of I(a;i) for 



■^Since a trivial loop gives an identity matrix, and taking the trace would give us a number 3. 
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Li and l{uj2) for L2, one can verify the second OPE in above formula. These wonderful OPE 
suggests that we got the right answer for the Wilson loops: the Wilson loop is classified 
by the irreducible representation of the gauge group, and the OPE is related to the tensor 
product of the corresponding representations. These facts play an important role in the 
physical approach to Geometric Langlands program in [35]. 

Let's now calculate the monodromy around the puncture, and the closed curve is shown 
in figure. 9. By taking a closed path representing the above closed curve in the monodromy 
graph, we have the following representation for the monodromy: 

Mp = Edge[e, 3].Face[y, 3].Edge[a, 3].Face[x, 3].Edge[6, 3].Face[y, 3]. 
Edge[c,3].Face[x,3].Edge[/,3].Face[y,3].Edge[(i,3].Face[x,3] (3.7) 

Notice that we use the different labels for the same edges {{e = c,a = f,b = d)), and at 
the end we should identify the coordinates of the same edge in the following way, say 

eo = ci,ei = cq. (3.8) 

we find a upper triangular matrices whose three diagonal terms (after making the matrix 
to have determinant one) are 

1 
aoaibobiCQCix'^y'^ 

Each diagonal term is the central element in the Poisson structure, and it is natural to 
regard the first two terms as the independent mass terms: 



711 = aoaibobicocixy, 722 = xy, 733 = ^ ^ ^ ^ _^ — 2~2- (^•^) 



mi = aoaibobiCQCixy , 1712 = xy. (3.10) 

We take the leading order term as the canonical map for this line operator with positive 
weights (as motivated by the rule in the Ai case): 

l{uj^) = aoaibobiCQCixy. (3.11) 

Here ujf = (1,0), and we use this rule following the study of the Wilson line. We can also 
take the least leading order term as the canonical map for the opposite label: 

aoaibobiCQCix^y^ 

here wj" = (—1,0). Similarly, by calculating the monodromy using the opposite orientation 
and take the leading order term, we have the canonical map 

I(wJ) = aoaibobiCQCix'^y'^, ^^2) = TT ' (3.13) 

aoaibobicocixy 

with UJ2 = (0,1), o;^ = (0,-1). The tropical a coordinates can be easily read from the 
canonical map. More generally, a closed curve around the puncture is labeled by a pair of 
integers {m,n), and the tropical a coordinates are 

Trop((m, n)) = Trop(wf "^'"^) + Trop(6J2^''^''^). (3.14) 

and the canonical map is just a single monomial in cluster X coordinates with the exponent 
determined by tropical a coordinates. 
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3.1.2 The missing line operator 

Completely same analysis can be applied to other simple closed curves, say M and C cycles. 
The conclusion is that the line operators supported on these closed curves are generated 
by the line operators in the fundamental and anti-fundamental representations, and they 
are labeled by the irreducible representations of SU(3) group. In the torus case, we have 
the generators Li, L2, Mi, M2, Ci, C2. If we take a duality frame where L cycle represents 
the weakly coupled gauge group, then Mj can be regarded as the 't Hooft loop operator. 

The above story seems to be in parallel with the Ai theory. However, something is 
missing here: first of all, all the line operators we found do not have the charges (tropical 
X coordinates) under the internal quiver nodes, and secondly there are only six generators 
for the loop coordinates and one should have at least 8 to match the total dimension of 
the cluster coordinates. It is natural that we could not find the extra ones by just looking 
at the closed curves since they should only carry the charge of the gauge group, and they 
are not able to probe the internal Coulomb branch operator carried by the T3 theory. 
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Figure 13: The OPE between the Wilson and t'Hooft line operators. 

Although a geometric picture for the missing ones is not available yet (we will turn to 
this later), we should be able to find them by doing the OPE of these usual Wilson and 
't Hooft line operators as they form a closed algebra. Let's do the OPE for Li and Mi 
(see the appendix for the canonical map of Mi), see figure. 13. In the OPE, it is natural 
to expect that the leading order component might be represented by the closed curve Wi 
shown in figure. 13, whose canonical map can be found using monodromy representation. 
We calculate the canonical map for Wi and subtract it from the product of I(Li) * I(Mi), 
and indeed we get a positive polynomial which means that it is indeed the top component 
of the OPE. We are left with a positive Laurent polynomial whose leading order term is 
(see appendix for the full Laurent polynomial): 



Kh 



a, 



2/3 1/3,1/3,2/3 2/3 1/3 



K'b 



-0 



-1 



xy + 



(3.15) 



This line operator is actually the missing one! The leading order monomial should give 
the tropical a coordinates for this new line operator. As usual. Let's find its tropical x 
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coordinates, and the answer is 



(0,0,0,0,0,0,-1,1) 



(3.16) 



so it indeed carries the charge under the internal quiver nodes, which strongly suggests that 
this is the missing line operator we are looking for. Similarly, one can find another internal 
line operator I2 by doing the OPE of L2 and M2 , and after subtracting the canonical map 
of W4, we find the canonical map 1(12) whose leading order term is 



^,j . 1/3 2/3,2/3,1/3 1/3 2/3 , 
I(/2) = Oq' a/ 6q' 6/ Cg' c^ xy + 



(3.17) 



and it has the tropical x coordinates 

(0,0,0,0,0,0,1,-1). (3.18) 

The coordinates of two new line operators are shown in figure. 14. With these two new line 



Tropical a coordinate 



Tropical x coordinate 







Figure 14: The tropical a and tropical x coordinates for /i and I2 which are found from OPE 
between the usual Wilson and 't Hooft line operators. 

operators, we have the following conjecture regarding the generators of the line operator: 
Conjecture: The line operators of A2 theory on torus with one full puncture are 
generated by the following eight nontrivial elements: 



Li,L2,Mi,M2,Ci,C2,Ii,l2; 



(3.19) 



By generators, I mean every nontrivial line operator with positive ^'^ tropical a coordinates 
can be formed by taking the product of them and take the leading order term in the OPE: 



I 



nGr% 



(3.20) 



where Gi are the above generators. The tropical data for these line operators are listed in 
table 1. 



^^If the initial tropical coordinates is negative, one can add the line operator with label (n, n) around the 
puncture to shift the a coordinates by a common positive constant, this will make the new a coordinates 
positive. So the non-trivial part is the line operator with positive tropical a coordinates, which is the same 
as A\ case. 
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Tropical a coordinates 


Tropical x coordinates 


Li 


a'^\'^h'^hY'^"'y"' 


(-1,0,1,0,0,0,0,0) 


L2 


aY'a^/h'JhY'x^lh'" 


(0,-1,0,1,0,0,0,0) 


Ml 


brhY'crcY'x^lh'" 


(0,0,-1,0,1,0,0,0) 


M2 


bY'bY'crcrxy^y^/' 


(0,0,0,-1,0,1,0,0) 


Ci 


a^/^a;/^cy^c^/\^/3,V3 


(1,0,0,0,0,-1,0,0) 


C2 


aY'aY'crcY'xy^y'/-^ 


(0,1,0,0,-1,0,0,0) 


h 


2/3 1/3,1/3,2/3 2/3 1/3 
«o «i K W Co c/ xy 


(0,0,0,0,0,0,-1,1) 


h 


1/3 2/3,2/3,1/3 1/3 2/3 
Oq' a/ 60' h^ Cg' c/ xy 


(0,0,0,0,0,0,1,-1) 



Table 1: The tropical a and a; coordinates for eight generators. 

3.1.3 Webs and Skein relations 

We found these extra line operators by doing OPE, which is purely algebraic as we use the 
canonical map. One may wonder if there is a geometric structure for them as line operators 
of 4d A/" = 2 theory should be represented by M2 branes wrapped on one skeleton on 
Riemann surface. We have already used all the closed curves, so new structure is needed. 
It seems the new structures are webs formed by three junctions, i.e. the Ii is formed 
by connecting two three junctions. Since the extra line operator has nontrivial tropical 
X coordinate on the internal quiver node, it is natural to put the junction at the internal 
quiver nodes, and it is useful to label the junction as black if the x coordinate is positive 
and white if the x coordinate is negative. Based on the tropical a coordinates of the line 
operator Ii and I2, it is then straightforward to find the a coordinate for the dual triangle 
of these two junctions, and the result is shown in figure. 15. 




Figure 15: Top: the junction and its tropical a coordinates on a dual triangle. Bottom: webs 
representing Ii and I2. 

Up to now, the junctions and webs are just a good way of representing the new extra 
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line operator. To make it as something really fundamental, we are going to do OPE for 
other line operators represented by the closed curves, and the result is shown in figure. 16. 
W2 (as well as W3) is represented by bipartite webs, and its tropical a coordinates and 
canonical map can be easily found: 

I(W2) = al/\y\hcl/'crx^/Y/' + ... (3.21) 
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Figure 16: Top: The geometric representation for the OPE between the Wilson and 't Hooft line 
operators. Bottom: the above OPE can be derived by using this simple Skein relation. 

Now a very strong check of the web structure is following: we can use another coordi- 
nate system by doing mutation on edge h and write down the tropical a coordinates of W2 
using the junction rule indicated in figure. 15, then we mutate back using the mutation 
sequences (/Li2,//4), (121,1^3), see figure. 17. The a coordinates of W2 in original coordinate 
system can be found using the mutation formula, and remarkably, the answer is the same 
as given by (3.21). One can check the transformation of the canonical map and they are 
actually related by the mutations! 

By examining the OPE shown in figure. 16, it is not hard to realize that they can 
be derived from a local Skein relation shown at the bottom of figure. 16. One might 
think that the OPE for any line operators can be found using this Skein relation. This 
is reasonable as one can do following: first draw the closed curves (or webs) for the line 
operators appearing in the operator product, and use the Skein relation to eliminate all 
four crossings, etc. However, the above story is not entirely correct, since there is a very 
crucial thing we want to address: there are trivial crossings between two closed curves, we'd 
better make sure such configurations are trivial and this would imposed the equivalence 
condition on the webs. By analyzing the equivalence conditions shown in figure. 18, we get 
the Skein relation in figure. 18, which we call as bubble reduction and square reduction 
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Figure 17: Top: the two quivers of two different triangulations are related by the quiver mutations 
{fi2,fJ'4:)i (a*1j/^3)- Bottom: The transformation of the tropical a coordinates of the W2 from one 
triangulation to another triangulation. 



respectively. So an irreducible web diagram has no two cycles and four cycles, and we 
conjecture that these are the full Skein relations. 



XX *'y—< 





o- 




Xk 



(-2) 




o 



Figure 18: A: The derivation of the bubble reduction and the square reduction. B: The full Skein 
relations of A2 theory. 

In the following, we will ignore the orientation of the edges in the web diagram, and it 
is easy to recover the orientations using the rules: the arrow is going from black junction to 
white junction. Using the above Skein relation, we have the following geometric picture for 
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line operators for A2 theory on once punctured torus: all the line operators can be realized 
as the bipartite webs without two cycles and four cycles on the torus^^. Let's look at two 
simple examples of using Skein relations to find OPE. 

Example 1: Let's consider the OPE of Ii and Li, and I2 and Li. Use the Skein 
relation, one have two line operators in the OPE, see figure. 19. 



I1*L1 



I2*L1 




L1*L2*M2 



Figure 19: Upper: OPE of I(/i) * I(Li). Middle: OPE of lih) * I(ii). Bottom: The OPE of 

I(Li)*I(L2)*I(A/2). 

Example 2: Let's study the OPE of L1L2M2 using Skein relation, and the result is 
shown in figure. 19: 



I(Li) * I(L2) * I(M2) = 1(A) + 1(B) + 1(C) + 3I(M2), 



(3.22) 



here A is the top component of 1(12) *I(Li), B is the top component of I(Li)*I(L2)*I(M2), 
C is the top component of I(L2) * I(Ci), and D is just M2. Notice that the coefficient 3 
before M2 is invisible from using the Skein relation. 

The above two examples are in perfect agreement with the calculation using the explicit 
canonical map. It is amazing that the simple Skein relation reproduces the correct OPE, 
which strongly suggest that the geometric picture we found for the general line operators 
are correct! 

3.1.4 Construction and reconstruction 

The web construction is based on the assumption that the line operators are built from eight 
generators (the leading component in OPE) . The construction of web from a coordinates 
is rather indirect: we need to first find the combinations of generators and then do OPE. 



^*Of course, if there is a marked point inside two or four cycles, they are allowed; we can only remove 
those contractible two and four cycles. 
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Moreover, we do not know what kind of constraints on a coordinates are needed to have 
a geometric web representation. So it is important to find a way of constructing webs 
directly from tropical a coordinate, which is also very useful for dealing with line operators 
on general Riemann surface. Indeed, we find similar reconstruction procedure as Ai case. 

Let's consider a triangle in the triangulation and assume the tropical a coordinates 
on the quiver nodes as ai; We further assume that the dual x coordinates of the internal 
quiver nodes to be positive. We multiple each tropical coordinates by 3, i.e. ao = 3ao, such 
that all the coordinates are integers now. 

The idea of reconstruction is to construct set of non-intersecting curves (junction) 
inside the triangle such that they give the desired a coordinates. There are actually seven 
kinds of curves: two curves around each vertex and one junction (only one type of junctions 
is allowed as we will see shortly). So let's assume there are ni (71,2) anti-clockwise (clockwise) 
curves around vertex A, and we use the similar notation for other two vertices; finally let's 
assume there are J black junctions inside the triangle. Summing up the contribution of each 
curve to the tropical a coordinate on seven quiver nodes, we have the following equations: 

2ni + n2 + mi + 2m2 + 2J = ao, ni + 2n2 + 2mi + m2 + J = ai, 

2mi + m2+h + 2/2 + 2J = 60, mi + 2^2 + 2li + h + J = h, 

2li + l2 + ni + 2n2 + 2J = cq, h + 2I2 + 2ni+n2 + J = ci, 

ni + 2n2 + mi + 2m2 + h + 2/2 + 3J = xq. (3.23) 

The solution is uniquely fixed by the tropical a coordinates: 



ni = -(ao + ci - xo), n2 = -(-ao + ai - 60 + xq), 








mi = -{ai + bo- Xo), m2 = oC-^o + h - cq + xq). 








h = (^1 + Co - Xo), h = ql-OO -C0+C1 + Xq), J = 


= ^(ao- 


- ai + 60 - 


-hi + co -ci). 

(3.24) 



Interestingly, J is the tropical x coordinates of the internal quiver nodes and therefore is a 
positive integer, which should be the case since only the junction can carry the tropical x 
charge of the internal quiver nodes. Notice that the number of marked points on one edge 
AB is equal to 

ni+n2 + J + mi+m2 = -{aQ + oi) = ai + oq, (3.25) 

O 

which actually only depends on the coordinate on the edge AB. Moreover, the number of 
curves coming into AB are ni + m2 + J = 2ao — ai , and the number of outgoing curves 
is n2 + m^i = 2ai — oq, which also only depend on the coordinates on the edge. To make 
all these numbers to be positive integer, we need to put the following constraints on the 
tropical a coordinates on the edge: 

ao + ai G z 

ao < 2ai, ai < 2oo. (3.26) 
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Moreover, we need the following constraints: 

a + /3 + 27Gz, 7<a + /3 (3.27) 

where 7 is the coordinate on internal quiver node, and a, /3 are quiver nodes in the same 
quiver triangle as 7, say (d, /3, 7) = (cq, ci, xq) in figure. 20. One can easily see that above 
three constraints would ensure the solutions in (3.24) are all positive integers! Similarly, if 
the tropical x coordinate of the internal quiver node is negative, we have the solution: 

'^i = o ("0 + ci - xo), 712 = -(co - ci - 61 + xo), 

^i = i,{ai + bQ- Xq), 7712 = -(ao - ai - ci + xq), 

h = o (^1 + Co - Xq), I2 = -{-ai + bo-bi + xo), J = -{-uq + ai - 69 + 61 - cq + ci), 

(3.28) 

and J is equal to the minus of the internal tropical x coordinate and therefore is positive. 
The constraints are the same as the above case. 

Example: Let's look at an example, the coordinates on the triangle are 

(ao, ai,bQ, 61, Co, ci,xq) = (5, 4, 1, 2, 3, 6, 5), (3.29) 

see figure. 20, so the marked points on each edge are {AB, BC, CA) = (3, 1, 3). Since the 
tropical x coordinates of internal quiver node is negative, we need to use solution (3.28). 
Substitute the coordinates into the formula, we have 

(m, 712, mi,m2,h,h, J) = (2, 0, 0, 0, 0, 0, 1), (3.30) 

so there is one white junction inside triangle and two clockwise curves around vertex A. 

Now given a set of a coordinates, we do the above reconstruction for each triangle, then 
we can form a set of closed non-intersecting webs and closed curves on Riemann surface. 
However, there are two complications we need to address: firstly the web has crossings 
if there are more than one three junctions inside the triangle; and secondly the marked 
points on the edge can be source (two arrows come out) or sinks (two arrows coming in). 
The first complication is easy to address: one replace each crossing by a piece with two 
junctions using the Skein relation, see an example in figure. 21. One always produce an 
irreducible sub web in this step as all the closed surfaces are hexagons. 

The second complication is a little bit difficult to deal with. First, let's color a marked 
points as black (white) if it is source (sink), and it is not hard to observe the following fact: 
the number of black and white vertices are the same. We can then connect the white and 
black vertices in pairs and produce a lot of crossings. The resolution of the crossings are 
different based on the origin of two crossing lines: 

• If the two lines coming from the vertices with different color, we keep the second term 
(two curves) in the Skein relation. 
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Figure 20: Top: construct curves from tropical a coordinates on a triangle. Bottom: construct 
webs on a torus from tropical a coordinates. 

• If the two lines coining from the vertices with same color, we keep the first term (two 
junctions) in the Skein relation. 

• If one of two crossing line is a ordinary line connecting two internal junctions, we 
keep the second term in the Skein relation; otherwise, we keep the second term. 

Using the above rules, we can produce an irreducible bipartite web around the edges. 
However, we do not know whether this procedure will produce an irreducible web in gluing 
the triangles together. This is true for all the examples we check, and it would be nice to 
have a proof though. 

In the above consideration, we start with a set of tropical a coordinates and then 
construct an irreducible web. On the other hand, if we start with a web, we would like to 
find its tropical a coordinates. This can be found by doing reverse Skein relation, namely 
one replace two junctions with a four crossings, see figure. 22^^. The above procedure 
stops if four legs of the two junctions are on the boundary. At the end, we are left with 
the basic building blocks: closed curves and Ii, l2- By taking a triangulation, one know the 
cluster coordinates for each individual component, and the coordinates are just the sum of 
that of components. 

For example, according to above method, the web diagram in figure. 22 is reduced 
to three generators (L2,L2,l2)- So the original web is the top component of the product 
I2*L2*L2 and its a coordinates are equal to Trop(l2)+Trop(L2)+Trop(L2), which actually 
equal to the a coordinates we used in last example. The result is in perfect agreement with 
the one shown in figure. 20, in which we construct the same web from same tropical a 



^^ Although the ordering of dissolving junctions does not matter, in practice it is better to start with the 
junction at the corner to avoid producing some weird diagram. 
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Figure 21: Top: resolve the crossing of the three junctions inside the triangle. Bottom: resolve 
the crossing around the edges. 

coordinates. Moreover, once the set of generators are found, one can find the canonical 
map of initial web by doing OPE. Now we have established a one-to-one correspondence 




Figure 22: A: replace two junctions with the crossing. B: We replace a web with a set of generators 
from which the a coordinates of the initial web is read. 

between the constrained tropical a coordinates (see [3.26, [3.27]) and irreducible bipartite 
webs. This description is exactly the same as the Ai case. 

It is worth to remark that bipartite network on bordered Riemann surface appears 
in many different physical contexts: 4d AA = 1 superconformal field theory [36, 37], BPS 
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quiver for 4d 7\A = 2 theory [32, 28], and scattering amplitude of TV = 4 SYM [38]. The 
web diagram we found here is of a different kind; the equivalence conditions are different: 
we do not have square move, but we have bubble and square reduction. 

3.1.5 A further check: Poisson brackets 

Let's now verify that the basic line operators (eight generators) are indeed the coordinates 
for the cluster variety associated with the moduli space of the local system on a torus 
with one full puncture. The most convincing check is to find the explicit Poisson brackets 
of these line operators, and show that there are indeed three commuting Hamiltonian as 
suggested by the underlying Hitchin integrable system. 

The Poisson bracket of line operators can be easily found due to the remarkably simple 
Poisson structure for the cluster X coordinates, which takes the following simple form: 

[Xj,Xj] = eijXiXj, (3.31) 

here €ij is the antisymmetric tensor read from the quiver. It is very easy to calculate the 
Poisson bracket for two loop operators using the canonical map which is just a positive 
Laurent polynomial in Xi. First, for two monomials A and B, the Poisson bracket is simply 

[A,B] = {J2ci\x'b)A*B, (3.32) 

i 

here a^ is the tropical a coordinates for the monomial A, and B is the tropical x coordinates 
for B. Using the above noted relation, we could easily calculate the Poisson structure 
between the eight generators using the explicit canonical map (See appendix B for their 
canonical maps.), and the Poisson brackets are indeed closed: 

[Li,Mi] = -h^Mi + h, [Li,M2] = -\liM2 + Ci 

[Li, Ci] = ^LiCi - /2, [LuC2] = ^LiC2 - Ml 

[Li, Ii] = ^Li/i - L2M1 + C2, [L1J2] = -\Lih + L2C1 - M2 (3.33) 

[L2, M2] = --L2M2 + h, [L2, Ml] = --L2M1 + C2 

[L2, C2] = '^L2C2 - h, [L2, Ci] = ^LsCi - M2 

[L2, h] = \L2I2 - L1M2 + Ci, [L2, h] = -\L2h + L1C2 - Ml (3.34) 

The other brackets are very similar 

[Ml, Ci] = -^MiCi + Li, [Ml, C2] = -^M2C2 + h 

[Ml, Ii] = -^MiLi + L1M2 - Ci, [M1J2] = \miI2 - M2C2 + L2 

[M2, Ci] = -^M2Ci + I2, [M2, C2] = -\M2C2 + L2 

[M2, h] = \M2I1 - MiCi + Li, [M2, h] = - \M2I2 + L2M1 - C2 (3.35) 
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and 

[CiJi] = -\cih + C2M2 - L2, [Ci,l2] = ^Ci/2 - C2L1 + Ml 

[C2, h] = IC2I1 - C1L2 + Ml, [C2, h] = -\C2I2 + CiMi - Li 

[/i, h] = L2M1C1 - L1M2C2. (3.36) 

There are actually two central elements which could be found from the above Poisson 
structure. One of them is 



/1/2 - L1M2C2 - L2M1C1 + C1C2 + M1M2 + LiL 



2 



111 , , 

= 3 + mi + m2 + mim2 H ! 1 , (3.37) 

mi ra2 mim2 

and mi,m2 is expressed as the Fock-Goncharov coordinates as 

777,1 = xy, 7772 = «ofli&o^iCoCixy. (3.38) 

Since ttii and 777,2 are the central elements of the Poisson bracket, and the above polynomial 
is also a central element. The other central element is much more complicated: 

L1L2M1M2C1C2 + Ll + Ll + Mf + Ml + Cf + Cl- If -II + 3/iCiLi 
+ 3/2C2L2 + 3/1L2M2 + 3/2^iMi + 3I1M2C2 + 3I2M1C1 - 3L1M2C2 - 3L2M1C1 

- 2CIL2M1 - 2CfLiM2 - 2CiLIMi - 2C2L1MI - 2C2LIM2 - 2C1L2MI 

- ClhLi - CII2M2 - CII2L2 - CfhMi - LIC2I2 - LlhM2 - LlCih- 

LlhMi - MlCih - MII1L2 - MIC2/1 - M| /2^i + llLiMi + C1IIM2 + C2IIL2 + C2IIM1 
+ /IL2M2 + Ci/|Li + C1C2L1L2 + C1C2M1M2 + L1L2M1M2 + /1/2M1M2 + I1I2L1L2 + C1C2/1/2 
C2LlMf + CiLJMl + CIlIM2 + CiLJMi + CIL2MI + CfliMf 

- C1C2I2L1M1 - C1C2/1L2M2 - C1I2L1L2M2 - C1/1L1M1M2 - C2I2L2M1M2 - L1L2C2I1M1 



6 111 6 

xy x^y^aQafftQ^fcQcf aoaibobiCQCi x^y^aoaibobicoci x^y^aoai^o^icoci 

6 

+ aofli^ofticoci + GxyaQaibobiCoCi + 6x'^y'^ aoaibobiCQCi 



xyaoaibobiCQCi 

+ x^y^aoaibobiCQCi + x'^y^aQa^bob^CQCi = 

1 1 7771 ^^-2 N 1 2 1 "^1 2 

12 + 6( |-rr7iH I-7772H 1 ) H o h 77727771 H k 1 h mim2. 

mi 7772 m2 7771 77727771 77727771 7772 

(3.39) 

Although the above formula seems to be very complicated, it is actually very easy to 
find using the canonical map. The key observation is the following: the central element 
is represented by the closed curves around the puncture and their dual x coordinate is 
zero! So if we can find the combinations of the simple generators such that their tropical 
X coordinates are vanishing, then its leading component is definitely a central element! 
We actually have only two independent choices /1/2 and L1L2M1M2C1C2, which give the 
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leading term in the above central element. The full expression is found by doing OPE of 
I(/i) * I(/2) and the OPE has the following form: 

I(/i) * I(/2) = Central terms + X] 11 '^i'" (^"^^^ 

i 

and the polynomial (Ii * -^2 — ^ Hi ^r') ^^ ^^^ central element. We also find the following 
three commuting Hamiltonians: 

Hi = Li* L2, 
H2 = Li + L2, 

H3 = LJCiMi + LIC2M2 + hh + M1M2 + C1C2 + L1L2 
- {I2C2L2 + M2L2I1 + hCiLi + M1L1/2). (3.41) 

In fact, one can choose any linear combinations of Li and L2 since they both commute 
with H^. The Poisson bracket for eight coordinates Li,Mi,Ci,Ii are closed, so they are 
good coordinates for the Poisson manifold. We have found two central elements and three 
commuting Hamiltonians, therefore it is an integrable system as it should be. The action 
variable can be found such that the symplectic structure has the canonical form 



OJ 



"^dHiAdei. (3.42) 



The detailed form for 9 coordinates will be presented elsewhere. These explicit calculations 
provide very strong evidence that we have the right answer for 4d AA = 2 line operators! 

3.2 General N 

3.2.1 Basic building blocks 

Let's consider line operators of 4d theory defined using general Ai\i^i (2, 0) theory on once 
punctured torus. The M = 2 theory in one duality frame is described by gauging the 

diagonal of two SU(A^) flavor symmetries of Tjv theory. The T^v theory has -^ 

Coulomb branch operators, and the SU(A^) gauge group has (A^ — 1) Coulomb branch 
operators; There are also {N — 1) mass parameters carried by the full puncture, so the 
rank of the charge lattice is 

R=(N- 1){N - 2) + 2(A^ - 1) + {N - 1) = N^ - 1. (3.43) 

The cluster coordinates are found by Fock and Goncharov: first take an ideal triangu- 

lation of the once punctured torus and then put '-^ more quiver nodes inside each 

triangle. The total number of quiver nodes is equal to the rank of the charge lattice of 4d 
theory, for more details, see [32, 28]. 

The ordinary Wilson lines are represented by the closed curves on the Riemann surface. 
The tropical o coordinates for Wilson loop in fundamental representation oji and anti- 
fundamental representation uon-i are found by calculating the monodromy M and M~ 
around the simple closed curve using the similar combinatorial method as Ai and A2 case, 
say appendix A for more details. 
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Once the canonical map (the positive Laurent polynomial) is found, its tropical a coor- 
dinates are read from the leading order term. There are other fundamental representations 
whose tropical a coordinates and the canonical map can be found by doing the OPEs of 
the two basic fundamentals. A useful fact is that the canonical map of the representation 
(n, 0, . . . , 0) and (0, . . . , 0, n) can be found by taking the trace of the monodromy matrices 
M" and M"". For example, we can do the product I((^i) * I(i^i), where wi is the funda- 
mental representation, and the same tensor product from the group theory tells us that 
the OPE has two components (2, 0, . . . , 0) and (0, 1, . . . , 0), since we know the canonical 
map for representation (2, 0, ... , 0), we can easily find the canonical map for the second 
fundamental representation and its a coordinates. Using OPE, one can find the tropical a 
coordinates for all fundamental representation oji. 
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Figure 23: The tropical a coordinates for the Wilson line in the ith fundamental representation: 
the quiver nodes in each straight line have the same indicated tropical a coordinates. 

Instead of using the OPE to find the tropical a coordinates of Wj, one can actually find 
the answer using the following requirement: the tropical x coordinates on the ith quiver 
node of a and h edge should be ±1 while all the other x coordinates are vanishing. The 
answer is: 
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(3.44) 



namely, the a coordinates take the maximal value at the ith line, and decreases along the 
left and right direction, see figure. 23. So it is natural to represent such line operators as 
the closed curve going through the ith line in the quiver diagram. 

After finding the tropical a coordinates and the canonical map of the fundamental 
representations, one can find all other representations by doing the OPE, which has the 
same support as the tensor product of the group theory. The conclusion is: the Wilson 
line is classified by the irreducible representation of SU(A^) group. 
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Again, these ordinary Wilson operators around three one cycles are not enough. The 
extra line operators can be found by doing the OPE between the Wilson line operators 
Li and 't Hooft operators Mj ■^^j and the extra line operator is found as the sub-leading 
component in the OPE. There are a total of {N — 1)^ pairs and the A^ — 1 pairs (Lj and 
M^^i) which will give the line operators in the C cycle, therefore we have (A^ — 1)(A^ — 2) 
number of extra line operators which exactly match the number of internal quiver nodes 
inside two triangles. The canonical map of the extra line operator can be found using the 
Poisson bracket: One calculate the Poisson bracket of Li and Mj, and the result would 
have two parts with different coefficients, and these two parts are the canonical map of two 
line operators in the OPE. 

These extra line operators can be represented geometrically using three junctions. By 
doing the explicit OPE, we find that there are - — ^—^ — ^^ type of three junctions which 
labeled by three positive integers satisfying the following conditions: 

i+j + k = N. (3.45) 

It is easy to see that the total number of solutions are - — ^-^ — ^ which equals to the 
number of internal quiver nodes inside one triangle. For each solution, there are two 
types of junctions such that the internal quiver nodes has positive (negative) tropical x 
coordinate. 

For the black junction, the arrows are coming out of the junction for all three legs 
while the arrows are coming into the junction for white junction, and we need to label the 
junction by (i,j, fc), see figure. 24. In the A^ = 3 case, there is only one solution for the 
above equation, therefore there is no need to specify the label for the legs. 

Let's now describe the tropical a coordinates on the quiver nodes of the dual triangle 
of black and white junctions. Again, the rules are found by looking the OPE of the line 
operators Lj and Mj. For the black junction, the tropical a coordinates of the internal 
quiver nodes and the third edge are found as following: first let's draw two lines Lj and 
and Mi in the i and j directions of the three junction; secondly, sum up the tropical a 
coordinates of the line operators Lj and Mi , finally we need to subtract appropriate integer 
number for the above tropical a on edge BC such that they are the tropical coordinates of 
the A;th fundamental representation (this can always be achieved), moreover, one need to 
subtract integer numbers from tropical a coordinates of some of the internal quiver nodes, 
the pattern is summarized in figure. 24. Similarly procedure can be done for the white 
junction. 

The dual tropical x coordinate of the junction would have non- vanishing value on only 
one internal quiver node. Using the above explicitly described tropical a coordinate, we 
find the position of the quiver node which the junction lives: 

f (a,6,c) = (z,j,fc). Black junction 
I {a,b,c) = {j,k,i), White junction. 



^®The convention of Li has maximal a coordinates on ith quiver node counting from left to right, while 
Mi is based on the counting from bottom to top, see figure. 23. 
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Figure 24: Top: the black and white junction with labels (i, j, A;), and the meaning of the label is 
that the tropical a coordinates is maximal on zth node counting in clockwise direction for the black 
junction, and anti-clockwise direction for white junction. Bottom: the rule of finding the tropical 
a coordinates from using fundamental representation tOi and uj: first sum up the a coordinates of 
LOi and ujj and then subtract 2 for the green nodes, and 1 on the red nodes. 

One new feature is that one can have M junction with labels ii,i2, ■ ■ ■ , im such that 

ii+i2 + ... + iM = N; (3.47) 

Such general junction can be formed using a sequence of three junctions by un-merging 
junction into a sequence of three junctions, so there is essentially nothing new. 

We oriented all the Wilson line in the same orientations as the first fundamental 
representation. For the following purpose, it is natural to choose a different rule for the 
orientation: the ith fundamental representation with I < i < ^ ^^s the same orientation as 
wi, otherwise it is oriented in the orientation as (^tv-i and the label is changed to i = N — i. 
If i = Y and is an integer, the corresponding representation is real, and we do not assign 
any orientation to it. With this convention, all the curves have the label with value less or 
equal than y. 

Accordingly, one need to change the orientation and the label for three junctions. The 
rules are actually very simple: because of the constraint i + j + k = N, there is at most 
one edge with label bigger than N/2, and if there is such type of edge, one simply change 
its orientation and change the label to N — i. The two representations have different 
advantages: the first one is good for constructing the complicated webs using the basic 
junctions, while the second one is good for the Skein relations we are going to study in 
the following. It is not hard to transform between these two representations. The basic 
junctions and its tropical a coordinates for ^3 theory are given in figure. 25, and we build 
some webs using the basic junctions in figure. 25. 
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Figure 25: Top: the basic junctions and its dual tropical a coordinates of Aj, theory. Bottom: 
some webs of A3 theory. 

3.2.2 Skein relations 

After describing the basic geometric building blocks for the webs, it is straightforward to 
develop the Skein relations and have a geometric representation for the OPE. In the mean 
time, any 4d line operator would be represented by the bipartite webs on the corresponding 
Riemann surface. For the purpose of Skein relations, it is useful to use second labeling for 
the three junctions, i.e. all the labels have value less than or equal ^. Two legs of two 
different junctions with same label and opposite orientations (including the legs without 
orientation) can be connected, and one can have multi-junction, which is different from the 
SU(3) case in which one can only have the three junction. 

The Skein relation between two oriented lines can be easily found by doing the OPE 
between the ordinary Wilson and t' Hoof operators in various fundamental representations, 
and the result is summarized in figure. 26. There are more Skein relations involving bub- 
bles, squares, and hexagons which are derived from obvious equivalence relations shown in 
figure. 26. A full Skein relations can be written down for A^ = 4 case, since there is actually 
just one type of junction with three legs labeled by (1, 1,2). Instead of labeling the edge 
with numbers, we use double line to represent the legs with label 2, and there is no need 
to specify its orientation since the corresponding representation is real. The two types of 
three junction are shown in figure. 27A, moreover, we have the four junctions which can 
be formed by combining the three junctions. The Skein relation for resolving the cross of 
two lines are shown in figure. 28. The details for the general cases are tedious and we leave 
it to the interested reader. 

Again, one can start with arbitrary tropical a coordinates and construct the webs 
using a similar reconstruction procedure as the one developed for A2 theory. One the 
other hand, one can start with an irreducible webs and find its tropical a coordinates. 
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Figure 26: Skein relations for general higher rank theory, here we assume i + j < y to draw the 
orientation for the new generated legs. 





Figure 27: The Skein relation for resolving the crossing for SU(4) theory. 

Our conjecture is that there is a one-to-one correspondence between the geometric and the 
algebraic description. 
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Figure 28: The full Skein relations of SU(4) theory. 

4. General bordered Riemann surface 

We discovered the basic ingredients such as line segments and three junctions with their 
coordinates on the dual triangle from the study of line operators on once punctured torus; 
we also find Skein relations which are very powerful in finding the OPE. The above ele- 
ments are completely local and do not depend on the global structure of Riemann surface, 
therefore they should still be valid on any bordered Riemann surface. We have the follow- 
ing conjecture: 4d line operator is built as an irreducible bipartite web on the Riemann 
surface, and the OPE can be found using Skein relations, etc. Other things like tropical 
a coordinates and canonical map are very similar. In the following, we are going to study 
the webs on several interesting bordered Riemann surfaces, and we use A2 theory as the 
main examples (the other cases are very similar). 

4.1 Ttv theory 

T/v theory is defined by a sphere with three full punctures, and is a strongly coupled SOFT 
with SU(N) flavor symmetry. There is no nontrivial closed cycle for three punctured 
sphere, which means that there is no weakly coupled gauge theory description, so there is 
no ordinary Wilson line. 

There are -^ Coulomb branch parameters and we should have non-trivial line 

operators represented by the webs using three junctions. The simplest one is built by 
two junctions: for each black junction with i + j + k = N, one naturally has a white 
junction on the back of the sphere to form a closed web, see figure. 29A. There are a total 
of {N — 1){N — 2) of such webs as one can easily count. In figure. 29A, we show the 
simplest web operators built from two junctions for T3 theory, which are denoted as Ai 
and A2- They can also be represented by the web on a plane graph in which the puncture 
is represented by a hole without any marked point, see figure. 29A. 
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Figure 29: A: the simplest webs on three punctured sphere. B: the OPE between above two webs 
are found using Skein relations, and there are two non-trivial webs in the OPE. 

In figure. 29B, we use Skein relations to derive the OPE of Ai and A2, and two 
nontrivial line operators appear which are denoted as Ji and J2. We can further do OPE 
between Ai and Ji, and find four more new line operators formed by eight three junctions, 
etc. Therefore, our conjecture is that the space of line operator is the same as the space of 
irreducible bipartite webs on three punctured sphere. 

Given a triangulation and the cluster coordinates [32], the tropical a coordinates of the 
web can be found using the basic rule for three junction. The canonical map for them is 
nontrivial since here we can not use the monodromy of the ordinary Wilson line operator 
and OPE. We do know the leading order term which can be read from the tropical a 
coordinates of the junctions, and the sub-leading order term is derived using the OPE: 



2/3 1/3,2/3,1/3 2/3 1/3 



2/3 1/3,2/3,1/3 2/3 1/3 



X . . . 



^\-^lj — "0 '^1 "0 "1 ^Q ^1 -^y ^ "'O "'I "0 "1 ^0 ^1 
J, . N 1/3 2/3,1/3,2/3 1/3 2/3 , 2/3 1/3,2/3,1/3 2/3 1/3 

/(A2) = Oq Oi V ^1 Co H xy + a^^' a/ b^' 6/ Cg' c/ y... 



(4.1) 



It would be interesting to find the full canonical map for these two basic line operators. We 
also draw the webs Ji and J2 in the triangulation, and a web which would appear in the 
OPE between Ji and Ai, see figure. 30. It is also not too hard to find that Ai,A2,Ji, J2 
form closed Poisson brackets even though we do not know the full canonical map. Because 
the leading component of Ai * A2 and Ji * J2 have trivial dual x coordinates, we find two 
constraints and two independent coordinates, which matches with the dimension of phase 
space of Hitchin integrable system. The explicit Poisson brackets and the Hamiltonian will 
be presented elsewhere. 
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Figure 30: Top: the cluster coordinates for T3 theory. Middle: Ji and J2 webs. Bottom: another 
web found in the OPE between Ji and Ai. 

4.2 Fourth punctured sphere 

Let's consider 4d theory defined by a sphere with four full punctures. This theory is 
described by a SU(A^) gauge group coupled with two T/v theories in one duality frame. The 
geometric picture of line operators for this theory is pretty similar: they are represented by 
irreducible bipartite webs. The OPE between two Wilson lines are derived using the Skein 
relations, see figure. 31. One can find the tropical a coordinates of these web operators by 
taking a triangulation, and the canonical map of various web operators can be found too. 
We leave the details to the interested reader. 

We would like to make a simple remark about the correlation function of 2d Toda field 
theory using our result (using the 2d/4d correspondence, the line operator on the Riemann 
surface will also be the observables of 2d Toda field theory.). Since the fourth punctured 
sphere can be derived by gluing two three punctured spheres, and the gluing only uses the 
line operator around the puncture, see figure. 32, and internal web operator is not used in 
gluing. In the 2d CFT interpretation, this means that the correlation function can be built 
from three point function, and the gluing is simple, which is consistent with 2d OPE. The 
nontrivial part of Toda field theory (counterpart of higher rank field theory) is that the 
three point function is complicated, as is manifested by the non-trivial line operators. The 
nice thing about the gluing picture is that once we understand the quantization of three 
punctured sphere, we should be able to understand the quantization of Toda field theory 
defined on any Riemann surface. 

4.3 Disc vi^ith marked points 

Let's consider a disc with marked points, which describes 4d Argyres-Douglas theory [28]. 
The special requirement is: one need to ensure that the tropical a coordinates of boundary 
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Figure 31: The OPE between line operators on fourth punetured sphere using the Skein relations. 




Figure 32: The gluing of two three punctured spheres is achieved by using the line operator around 
the puncture. 

quiver nodes to be zero, as we only counter the internal quiver nodes in studying of 4d 
A/" = 2 theory. Another important rule is that for the curve connected with two boundary 
points, the canonical map is defined as 



1(0 



1 0) *M* I 
0, 



(4.2) 



here M is the monodromy around the path and we set the boundary coordinates to be 
one in the end, so the answer does not depend on the boundary coordinates. 



■^^There are choices on which edge we should start or end in the monodromy graph, here we start and 
end with the edge which have the same orientations as the curve. 
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For simplicity, let's consider A2 theory on a disc with four marked points, which 
actually represents D4 Argyres-Douglas theory [28]. A lamination satisfying the above 
constraint is shown in figure. 33, and its tropical a coordinates is also shown. The canonical 
map is very simple I(Ti) = XiX^, which is actually the central element in the Poisson 
structure (consider only the quiver nodes insider the quadrilateral). Another lamination 
using web is shown in figure. 33, and the canonical map of W and the other two curves 
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SO the canonical map of this lamination is 

1{T2) = 1{W) * I{A) * I{B) 



1 1 1 

Xx X1X4 X1X3X4 



(4.3) 



(4.4) 







Figure 33: Left: two laminations on fourth punctured disc of A2 theory. Right: The tropical a 
coordinates of these two laminations. 

The same theory can be constructed using Ai theory on a disc with four marked 
points and one bulk puncture, and the triangulation which give the same quiver as A2 
triangulations are shown in figure. 34. The lamination with the same tropical a coordinates 
of T2 is found using the reconstruction procedure described in section 2, see figure. 34. 
The canonical map of this lamination T3 is found by using the monodromy representation 
of the SU(2) local system described in section 2: 

_i _i _i _i _i 111 11 _i 1 _i _i 

2 



1(A) = X,\ 1(B) = X, 'X, 'X, 'X, ^ 1(C) = X^XiXi+X.'XiX, '+XiX, 'X, 



(4.5) 



^*One first find the canonical map of curve connecting two opposite boundaries using the monodromy, 
and then do the OPE to find canonical map of web operator. 
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and the canonical map for this lamination is 
I(r3) = 1(A) * 1(B) * 1(C) = 



1 1 1 

Xi X1X4 X1X3X4 



(4.6) 



Remarkably, this is equal to the answer found using SU(3) representation. This result is 
clearly a very strong evidence that our generalization of lamination to higher rank theory 
is correct. 



T3 



















x^ 


/" 


\ 




,y^ 


\; 


y 




Xl^^ 


/-..x 




Figure 34: A A\ lamination which gives the same tropical a coordinates as the lamination T^. 



5. Conclusion 

In this paper, we initiated a systematical study of half-BPS line operators of 4d AA = 2 
theory (theory of class 5) defined using full defects. Although the basic building blocks 
like line segments, three junctions are described for any A^, we mainly focused on A2 and 
Aj, theories, it would be interesting to further study general cases, i.e. giving the full Skein 
relations, finding the construction and reconstruction procedure as we did for A2 case, etc. 
There are many other open questions and further generalizations: 

• We found the canonical map of line operators using OPE, which is straightforward 
but tedious. It would be very useful if the answer can be directly found using web 
picture. This is important for T/v theory as the OPE method is not applicable. 

• It would be interesting to further study the OPE and general properties of the in- 
teger coefficients like its support and value. These results are useful in checking the 
conjectures given in [39]. It is also interesting to compare our results with the OPE 
of Wilson-'t Hooft operators studied in the context of 7\A = 4 SYM [35, 40, 41, 42] 
and three dimensional theory [43]. 

• The non-full puncture case is worth to pursue, since in this case there are theories with 
Lagrangian descriptions and many field theory results are available [44, 45, 46, 47] 
(see also some results using Toda field theory [48, 49]). The basic tools like junctions. 
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Skein relations should still be valid, and the cluster coordinates are given in [32, 28] 
which can be used to parameterize the line operators, etc; The main difficulty is that 
there are no monodromy calculation tools available, so we don't know how to get the 
canonical map, currently we are trying to deal with this problem 



• Finally, it would be interesting to study various physical quantities related to the line 
operators, i.e. Superconformal index (see [51] for some Ai cases), framed BPS states, 
S"^ partition function in the presence of line operators (see [52, 53] for Ai theory 
with Lagrangian description) etc; it is also interesting to study how the line operator 
labels change under the pants decomposition, which would be a nice justification of 
S duality. 

One actually need a quantum version of line operators to answer many questions raised 
above. The quantum deformations of the line operator can be studied using quantum 
cluster algebra [39], for the Ai case, see [27, 10, 54]. A complete treatment about these 
quantum line operators of higher rank case will be given in a future publication [55] . 

In this paper, we mainly focus on A lamination and do not touch X lamination. The 
main difference is the webs of X lamination can end on bulk punctures. Such X lamination 
is important in representing 4d BPS states in the Ai case [56]. Naively, we expect the higher 
rank X lamination should also give the 4d BPS states, and such interpretation matches 
perfectly with the expectation from the BPS quiver approach in some special cases we 
check, and we hope to come to this point in the near future. 

The three junctions play a crucial role in our construction, and our derivation is based 
on the canonical map, cluster coordinates, etc. It would be interesting to see whether it 
can be realized as the one quarter BPS surface operator of (2, 0) theory. Similar web-like 
objects in M theory have been studied in [57, 58, 59], it is interesting to see if there is any 
relations between the junctions (superficially it seems that the number of junctions in our 
case is smaller than theirs.), see also the quarter BPS supergravity solutions related to M2 
branes in [60]. If these webs can indeed be realized as quarter BPS surface operators, we 
expect them to be important in understanding the dynamics of (2, 0) theory. 

These line operators would have a lot of applications in various physical contexts due 
to the magical six dimensional construction, see figure. 35. It seems that all of them can 
be pursued using our result, and it would be exciting to do some concrete calculations. 

Our work also has a lot of direct relations to mathematics. Firstly our web structure 
seems to be the geometric picture for the higher lamination proposed in [27, 39] (see also 
recent work on geometric realization of higher lamination using affine Grassmannia [61]); 
The web diagrams studied by us are also called tensor diagrams which are used to study 
cluster algebra in [62]. Secondly it is amusing to note that our Skein relations coincide 
with the one used in the study of Khovanov homology, see [63, 64], it would be interesting 
to see if there is any deep connection. Thirdly the OPE of Wilson and t'Hooft loops [35] 
and quantization of Hitchin integrable system [65] play an important role in the study of 
Geometric Langlands duality. Our results are both related to the OPE of line operators 
and quantization of Hitchin integrable system, and it is interesting to see if there are any 
potential applications. 
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Figure 35: The use of line (web) operators. 
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A. Calculating monodromy 

The cluster coordinates for higher rank theory defined using full punctures are found in 
[27]. It was conjectured that they are the BPS quiver for corresponding A/" = 2 theory 
in [32] which also gave a brane construction. One start with an ideal triangulation of 
bordered Riemann surface, and then further triangulate triangles in the triangulation. 

There are -^ new quiver nodes inside the triangle with three vertices ABC, and 

each internal quiver node can be labeled by three positive integers satisfying the condition 

a + b + c = N, (A.l) 

here a is defined as the distance to the dual side of the vertices A, etc, see figure. 36. 

Let's now summarize the ingredients for calculating the monodromy of higher rank 
local system following Fock and Goncharov [27] . For each triangulation, one can construct 
a dual graph: put one trivalent vertex inside each triangle and connect those trivalent 
vertices together to form a dual graph. We form a hexagon around each vertex and a 
square around each edge in the dual graph. We call dual graph together with the hexagons 
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S: Cross edge of dual graph 



F: Inside triangle 



E; Parallel with edge of dual graph 



Figure 36: Top: the monodromy graph. Bottom: the label for coordinates inside a triangle, here 
the first index means the distance to the dual edge of vertex A. 

and squares as the monodromy graph. There are three types of edges in monodromy graph: 
S edge which intersects with the dual graph, E edge which is parallel with the edge of dual 
graph, and F edge which is inside the triangle, see figure. 36. We choose an orientation for 
edges in monodromy graph: the orientation is taken to be anti-clockwise for edges circling 
around the punctures, and all the edges in hexagon is taken to be in clockwise orientation. 
The monodromy is calculated by assigning a N x N matrix to each edge in monodromy 
graph. The E type is a diagonal matrix depending on the coordinates on the edge of the 
triangulation with which edge E intersects: 



E{a)[N] = diag(aoaia2...aAr-i, ai02...aAr_i, . . . , a^-i, 1) 



(A.2) 



where the starting point for the label is taken as the vertex where edge E circles around, 
see figure. 37. Notice that the matrix for other edge parallel with E is denoted as E{a) 
since we have different labeling for the coordinates. The coordinates a are related to the a 
coordinates as 

(A.3) 



ao = a7v_i, , 



• ^"-i 



aN~j~ii ■ ■ ■ ,^N~i 



ao. 



The S Matrix also has a simple form: 



(-l)i-i j = N + l-i 







(A.4) 



The F(ace) matrix is more complicated and it is composed by Fi and Hi matrix. The Hi 
matrix is a diagonal one which has the form depending on position of the internal quiver 
node: 

Hi{x) = diag{x,x,x,l,. . .,1), (A. 5) 
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and Fi matrix is defined as 



(Fi. 




(A.6) 



So Fi matrix is a lower triangular matrix with entry one on the diagonal. Using the above 
two building blocks, the F(ace) matrix is expressed as 



1 N-2 

F[x,N]= J] { H {H,+i{XN-i-2,i-jj-i)F^)FN- 

j=N-l i=j-l 



(A.7) 



Notice that here an internal quiver node is labeled by three integers satisfying a + b + c = 
N — 3, which is given by subtracting one from our previous labeling. For small A^, the 
matrices are 



F[x, 2] = Fi, F[x, 3] = F2H2iXo^o,o)FiF2, 

F[x,4] = F3H3{Xo^i,o)F2F3H2{Xi^o,o)FiH3{Xo,o,i)F2F3. 



(Ai 



There is an easy way of remembering this sequence by looking at the quiver diagram in 
figure. 36. If the monodromy path is going from edge AB to AC, then the above product 
is decomposed into A^ — 1 steps. In each step, each quiver node on the line parallel to edge 
AC contribute a factor Hi+iFi, where i is the distance to vertex A of that quiver node. In 
each step, the sequence of product is to start from top to bottom, see figure. 37 for the 
arrows indicating the order of product; finally one need to add a matrix Fat-i at the end 
in each step. 





Figure 37: Top: the labeling for the coordinates on the edge with which E intersects, and this 
labeling is used in the calculation of the monodromy. Bottom: the rule of remembering the formula 
for calculating the F matrix going from AB to AC. 

One should be really careful about the fact that the labeling of quiver nodes are quite 
different if the path is going through the other vertex of the triangle, so it is better to use 
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different letters for the face matrix inside the same triangle, say X,X,X, and finally one 
should relate three set of coordinates by looking at their positions, say Xa,b,c = Xb^a,c, etc. 
The above matrices are the basic ingredients for calculating the monodromy, and if a 
path goes in opposite direction as the chosen orientation, we just use the inverse of above 
matrices. The monodromy around the square and the hexagon is an identity matrix up 
to a constant scalar, therefore the monodromy only depends on the homotopy class up 
to a complex scalar, and it defines a PGL(A^, C) local system. The monodromy for the 
SL(iV, C) local system is defined as 

/(/) = |Det(M)-i/^Tr(M)|. (A.9) 

here we take the absolute value to make sure that all the expressions are positive. 

B. Collection of canonical maps 

For A2 theory defined on once punctured torus, there are eight nontrivial generators. The 
quivers and the geometric representations for these eight line operators are shown in fig- 
ure. 11 and 15. From the quiver, one can read the Poisson brackets for cluster coordinates, 
and one then use them to calculate the Poisson structure between the line operators com- 
bining the explicit expression of canonical map. We collect canonical maps for these eight 
generators, and the interested reader can do all kinds of calculations presented in the paper. 

1 b'/' X1/36V3 xV3„l/3ftl/3 
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